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Abstract. We present a microscopic model for neutrino induced one-pion production off the nucleon and its implementation 
for the purpose of calculating coherent pion production in nuclei. We further criticize the use of the Rein-Sehgal model for 
coherent pion production by low energy neutrinos. In particular, we show how the approximations in that model give rise to a 
much flatter differential cross section in the r/ = E n (l — cos Q n ) variable. We discuss the limitations intrinsic to any approach 
based on the partial conservation of the axial current hypothesis and the inability of such models to properly determine the 
angular distribution of the outgoing pion with respect to the direction of the incoming neutrino. We show the effects of those 
limitation for the case of the ^ differential cross section. 

Keywords: Neutrino reactions, coherent pion production 
PACS: 25.30.Pt, 12.15.-y, 13.15.+g 

INTRODUCTION 

Neutrino induced one-pion production off nucleons and nuclei in the intermediate energy region is a source of relevant 
data on hadronic structure. Pions are mainly produced through resonance excitation and these reactions can be used 
to extract information on nucleon-to-resonance axial transition form factors. Besides, a proper understanding of these 
processes is very important in the analysis of neutrino oscillation experiments since pion production could be a source 
of background in those experiments [1,2]. 

In reactions on nuclei, pions can be produced incoherently or coherently. In the latter case the nucleus remains in its 
ground state. Coherent reactions are controlled by the nucleus form factor and are more forward peaked than incoherent 
ones. Experimental analyses of the coherent reaction rely on the Rein-Sehgal (RS) model [3] which assumes that 
coherent pion production is dominated by the divergence of the axial current and can thus be related to the pion- 
nucleus coherent scattering through the partial conservation of the axial current (PCAC) hypothesis. For instance for 
coherent 7C° production the RS model approximates the coherent cross section for both neutrino and antineutrino 
induced processes by 
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with G the Fermi constant, M the nucleon mass, E the incident neutrino energy and f n the pion decay constant. Besides 
q is the lepton four-momentum transfer, x = —q 2 /2Mq°, y = q°/E, and t — (q — k n ) 2 with k n the pion four-momentum. 
t equals — (q — k n ) 2 if, as usually assumed, the nucleus recoil is neglected (q° = k° % = E n ). F^(t) is the the nuclear 

form factor calculated as F^(t) — f d^r e x ^ klt ^' r {p p {r ) + p n (r )} with p p ( n ) the nuclear proton (neutron) density, 
normalized to the number of protons (neutrons). f a b s is a f— independent attenuation factor that takes into account 

the distortion of the final pion. Finally da ^ n Q is the differential pion-nucleon cross section evaluated 

at t = 0. The t = approximation in the pion-nucleon cross section is not needed and it can be justified only if the 
nuclear form factor is sufficiently forward peaked. The larger the pion energy and the heavier the nucleus, the better 
this approximation becomes. In the original paper [3], the model was applied to medium size nucleus, aluminum, and 



neutrino energies above 2 GeV, for which the relevant pion energies are quite high. However, as we pointed out in [4], 
for neutrino energies below 1 GeV and lighter nuclei, like carbon or oxygen, the nuclear form factor is not enough 
forward peaked to render the finite t— dependence of the pion-nucleon cross section negligible. The distortion factor 
F a b s is also an oversimplification since in any realistic scattering model this factor should depend on t . The recent work 
of Berger and Sehgal [5] already corrects some of the problems in the RS model. 

In Ref. [6] we discuss in detail the approximations inherent to any PCAC based model and, in particular, the 
approximations in the RS model. As we show below, the neglect of non-PCAC terms and the implicit assumption 
(only correct for q 2 — 0) by PCAC based models that g^g^- — = Jk ch%dt > wriere tyk n q is the azimuthal pion angle 
in a plane perpendicular to q, affects the determination of the differential cross sections with respect to the angle 
made by the pion and the incident neutrino, leading to flatter distributions in that variable. In the particular case of 
the RS model the t — approximation enhances this unwanted effect. This could have implications for the recent 
determination by the MiniBooNE Collaboration of the neutral current (NC) coherent K° production rate [7]. 

Our approach to the problem, as others in the literature, does not rely on PCAC. We use a microscopic model in 
which coherent pions are mainly produced by virtual A-hole excitations in the nucleus, as well as additional processes 
that are required by chiral symmetry. The model includes the modifications of the A propagator in the nucleus by 
means of a sophisticated evaluation of the A self-energy in the medium, and it treats the final pion distortion in a 
realistic way by solving the Klein-Gordon (KG) equation for a pion-nucleus optical potential. 



MICROSCOPIC MODEL AT THE NUCLEON LEVEL 



The model that we use to describe the elementary production process at the nucleon level is the one derived in Ref. [8]. 
In Fig. 1 we give the different contributions for the charged current (CC) W + N —>N'n case (details for the antineutrino 
induced process and the NC case can be found in [8]). In addition to the dominant A pole (AP) mechanism (weak 
excitation of the A(1232) resonance and its subsequent decay into Nn), the model also includes the crossed A pole 
term (CAP) and other background terms that are required by chiral symmetry: Direct and crossed nucleon (second 
row) pole terms ( NP, CNP), contact (CT) and pion pole (PP) contribution (third row) and finally the pion-in-flight 
(PF) term. In Ref. [8], we found that background terms produced significant effects in all channels. The least known 
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FIGURE 1. Model of Ref. [8] for the W + N — > N'n reaction. The circle in the diagrams stands for the weak transition vertex. 

ingredients of the model are the axial nucleon-to-A transition form factors, of which Cj gives the largest contribution 
(See Eq.(l) of Ref. [9] for a form factor decomposition of the NA weak current). Besides, within the Adler model that 
we use [10], Cj determines all other axial form factors. This strongly suggested the readjustment of that form factor 
to the experimental data, which we did by fitting the flux-averaged v^p — » fJ,~pn + ANL [11, 12] q 2 -differential cross 
section for pion-nucleon invariant masses W < 1.4 GeV, for which the model should be appropriate. Assuming the 



shape [13] 
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we found C% (0) =0.867 ±0.075 and M AA = 0.985 ±0.082 GeV. Our full model, thus obtained, lead to an overall better 
description of the data for different CC and NC, neutrino- and antineutrino-induced, one-pion production reactions off 
the nucleon [8]. The Cj (0) value is significantly smaller than the traditionally used value of about 1.2 deduced from 
the off-diagonal Goldberger-Treiman relation. This reduction of Cg (0) is consistent with recent results in lattice QCD 
[14], quark model [15] and phenomenological studies [16]. 



COHERENT PRODUCTION IN NUCLEI: MODEL AND RESULTS 

Now the process consists of a weak pion production followed by the strong distortion of the pion in its way out of the 
nucleus. In the coherent production the nucleus is left in its ground state and the hadronic amplitude is written as a 
sum over the amplitudes for each nucleon. Neglecting for the moment nonlocalities and pion distortion one can write 
for the hadronic amplitude 

- / d*re^- k *)- T { Pp {r)[j?£ n {q,k n )\ +p n {r) [fUl,**)] } O) 

where the approximation q° = k () K (neglect of nucleus recoil) has been used. One can see from Eq.(3) that the coherent 
pion production process is sensitive to the Fourier transform of the nuclear density for momentum q — k n . J?^k il^x) 
in Eq.(3), stands for the nucleon spin averaged WN — > N7Z or Z°N — > Ntz weak transition amplitude 

f£ K (q,k K ) = Wtir(p')r£ NK ii r (p) i = AP, CAP, NP, CNP, CT, PP, PF (4) 
i r p 

where the u's are Dirac spinors for the nucleons, normalized such that uu = 2M, and the four-vector matrices T^. Nn 

can be read from the explicit expressions of the pion production amplitudes given in Ref. [8]. p and p' — p + q — k K 
are the initial and final three momenta of the nucleon. Those momenta are not well defined and we approximate the 

four- momentum of the initial nucleon by p^ = \J M 2 + \(k n — q) 2 , )_ Hence we assume that the initial nucleon 

momentum is (k K — q)/2 and the final one is — (k K — q)/2, with both nucleons being on-shell. The momentum transfer 
is equally shared between the initial and final nucleon momenta. Setting p = —p' = (k K —q)/2 eliminates some non- 
local contributions, and it greatly simplifies the sum over all nucleons, which can be cast in terms of the neutron 
and proton densities. Furthermore, the sum over spins in Eq. (4) can be also easily performed for p = —p' since 
u r{p' = —p) = y°u r (p), so that 



^u r (p' = -p)T» m Ur{p) = iTr((^+M)y°r^J, i = AP,CAP,NP. . . 



(5) 



Given the importance of the A— pole contribution and since the A properties are strongly modified inside the nuclear 
medium we consider some additional nuclear corrections to this contribution to include the effect of the self-energy 
of the A in the medium L^(p(r)). We follow the same approach as in Ref. [17], which is based on the findings of 
Refs. [18, 19, 20] . Thus in the A-propagator, we make the substitutions M A ^M A + Re£ A and T A /2 -> rP au!i /2 - 
ImEA and take Ea(p (r )) and r^ aull /2 as explained in Sect. II-B of Ref. [17]. Once these corrections are included 
depends on r as well and the nuclear form factor can no longer be factorized out as in the RS model. 

So far the formalism has used the bound wave functions of the nucleons in the nucleus, which appear via the proton 
and neutron densities, and has considered only a plane wave for the pion. Pion distortion effects are important, specially 
for \k n \ < 0.5 GeV, and are considered by replacing in Eq. (3) 

e-^^yZinh) (6) 
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FIGURE 2. Effects of A in medium, % distortion and nonlocalities in % momentum on the 4& differential cross section evaluated 

ClK 7Z 

for the v 16 -» ^ l6 Olt+ reaction at E v = 0.6GeV. 



The pion wave function (p*(r;k K ) corresponds to an incoming solution of the KG equation, 



-v 



2E n V opt (r) <?* K (r,k n ) = E^(r;k K ) , 



(8) 



with V p t (r) the optical potential which describes the ft-nucleus interaction. This potential has been developed micro- 
scopically and it is explained in detail in Refs. [19, 20]. This theoretical potential reproduces fairly well the data of 
pionic atoms (binding energies and strong absorption widths) [19] and low energy ^-nucleus scattering [20]. At low 
pion energies, it is an improvement over the one used in [17, 21], that was based on A dominance of the %N interac- 
tion. The replacement in Eq. (7), that takes into account the fact that the pion three momentum is only well defined 
asymptotically when the pion-nucleus potential vanishes, induces some non-localities in the amplitudes. To treat these 
non-localities we have adopted the following scheme: 

• In the AP, CAP, NP, CNP terms, we note that there exist either a NNn or a NAtz vertex (see Eq. (5 1 ) of Ref. [8]), 
which induces a factor k" in the amplitudes. Indeed, for those terms we could re-write 



fg Nn {r,qM = {k K ) a f^{r,q,k K ), i = AP,CAP ,NP ,CNP . 



(9) 



We do not consider any non-locality in the tensor 



and we use the prescription of Eqs. (6) and (7) to account 



for k n in the contraction between k" and JffjSL in Eq. (9). For k ( l we shall use the asymptotic pion energy. This 
approach to treat the non-localities is equivalent to that assumed in Refs. [17, 21]. 
• We do not consider any non-locality for the CT and PP contributions. 

The effects of the modifications of the A properties in the medium, pion distortion and nonlocalities in the pion 
momentum are shown in Fig. 2. As can be seen all corrections are very relevant. Very recently, it has been claimed that 
nonlocalities in the nucleon momentum could be as important, giving rise to a large reduction of the coherent cross 
section [22]. The calculation in Ref. [22] has been done without A in medium corrections or pion distortion and further 
studies are needed to clarify this point. 

As mentioned before, background terms gave an important contribution at the nucleon level [8] but they turned out 
to be irrelevant for coherent production in symmetric nuclei [4], and thus the reduction we found in C5 (0) amounts 
to an important decrease in the coherent cross section. In Fig. 4 we give cross sections for CC and NC coherent pion 
production on carbon and oxygen. Other observables are discussed in Ref. [4]. 

Let us here concentrate on the ^ differential cross section with rj = £^(1 — cos the variable introduced by the 

MiniBooNE Collaboration in its analysis of NC coherent 7T° production [7]. In Fig. 4 we show results for NC coherent 
n° production on carbon at E v = 1 GeV. In the left panel the results have been evaluated without distortion or any 
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FIGURE 3. Muon neutrino/antineutrino CC and NC coherent pion production off nuclei from carbon (left) and oxygen (right) 
targets as a function of the neutrino/antineutrino energy. 
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FIGURE 4. Comparison of different predictions for the ^ differential cross section for NC coherent K° production on carbon 
at E v = 1 GeV. See text for details. 



in medium correction. We compare the results obtained with the RS model, the one used by experimentalists in their 
analyses, with the results obtained in the microscopic model of Refs.[8, 4]. In the latter case and for simplicity, we 
have only taken the dominant AP term and, in order to fix the normalization, we use Cj (0) = 1.2 which is the value 
implicit in any PCAC-based model. Besides we show, with triangles, results of the microscopic model eliminating all 
non-PCAC contributions. Comparison of the latter with the RS model shows the effect of the t = approximation 
used in that model. Cross sections are larger in the RS model and what is more important this model produces a 
much flatter distribution (see line "Rein-Sehgal rescaled"). As described in Ref. [6], one can easily eliminate the t = 
approximation in the RS model by using the full non spin-flip part of the pion nucleon elastic cross section. The results 
for this case are shown under "Modified Rein-Sehgal" with a dashed line. We see the agreement with the only-PCAC- 
terms microscopic model is very good in this case. Another important piece of information can be extracted from 



the comparison between the full microscopic calculation and the only-PCAC-terms one. We see the latter produces 
a flatter distribution. This is mainly due to the 21cdy3tdf = 27? dxdydt approximation inherent to any PCAC-based 
model. In the right panel we show results that include distortion and in medium corrections. Again we see the RS 
model produces larger and much flatter distributions than the microscopic calculation due to the t = approximation 
and the inadequacy of the final pion distortion used. Besides the effect of the dxdyd f dt j, k — dx d dydt approximation 
can also be seen by comparing the two microscopic calculations. 



CONCLUSIONS 

To conclude we would like to stress the following: First, background terms are relevant at the nucleon level and 
should be included in any microscopic calculation. Even though they play a minor role in coherent production, their 
inclusion affects the determination of the C5 (0) value which is a crucial quantity for coherent production as it fixes 
the strength of the dominant AP contribution. Second, the Rein-Sehgal model was cleverly devised for its use at high 
neutrino energies and medium/large nuclei but we think we have shown it is inadequate for low neutrino energies 
and small nuclei. The use of the RS model in experimental analyses could affect the determination of the rate of 
coherent production. We believe this is the case for the recent determination of NC ft coherent production by the 
MiniBooNE Collaboration [7]. Last, some of the problems of the RS model are shared by all PCAC-based models: 
the neglect of non-PCAC terms and the implicit assumption dxdy d dtd<Sfk = 5^ dxdydt ma kes those models unable to 
properly determine differential cross sections with respect to any variable that depends on the angle made by the pion 
and the incident neutrino. We have shown that for the case of the differential cross section. 
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